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Abstract 

We have substantially extended the high-temperature and low-magnetic-field (and the related 
low-temperature and high-magnetic-field) bivariate expansions of the free energy for the conven- 
tional three-dimensional Ising model and for a variety of other spin systems generally assumed to 
belong to the same critical universality class. In particular, we have also derived the analogous 
expansions for the Ising models with spin s = 1,3/2,.. and for the lattice euclidean scalar field 
theory with quartic self-interaction, on the simple cubic and the body-centered cubic lattices. Our 
bivariate high-temperature expansions, which extend through i^^^, enable us to compute, through 
the same order, all higher derivatives of the free energy with respect to the field, namely all higher 
susceptibilities. These data make more accurate checks possible, in critical conditions, both of the 
scaling and the universality properties with respect to the lattice and the interaction structure and 
also help to improve an approximate parametric representation of the critical equation of state for 
the three-dimensional Ising model universality class. 

PACS numbers: 05.50. +q, 64.60.De, 75.10.Hk, 64.70.F-, 64.10. +h 
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I. INTRODUCTION 



We present a brief analysis of high-temperature (HT) and low-field expansions for the free 
energy of the conventional 3D Ising model in an external uniform magnetic field, extended from 
the presently available^."— order 17 up to order 24 in the case of the simple-cubic (sc) lattice, and 
from the order 13 up to 24 in the case of the body-centered-cubic (bcc) lattice. In addition to 
the conventional Ising model (i.e. with spin s = 1/2), we have considered also a few models with 
spin s > 1/2, and the lattice scalar euclidean field theories with even polynomial self-interaction. 
All results for the simple Ising system in a field can be readily transcribed into the lattice-gas 
model language and therefore are of immediate relevance also for the theory of the liquid-gas 
transitioii^i^. The HT and low-field expansions of the spin-s Ising models can be transformed^i^^ 
into low-temperature (LT) and high-field expansions. 

The spin-s Ising model in an external magnetic field H is described by the Hamiltoniaitii"— 

<ij> i 

where Si = —s,—s + 1, ...,s is the spin variable at the lattice site i, m is the magnetic moment of 
a spin, J is the exchange coupling. The first sum extends over all distinct nearest-neighbor pairs 
of sites, the second sum over all lattice sites. The conventional Ising model is recovered by setting 
s = 1/2. 

The one-component self-interacting lattice scalar field theory is described by the 
Hamiltoniani^"^— 

n{<P} = -Yl ^i^3 + Y.^y{^^) + H(t>^). (2) 

<ij> i 

Here — oo < (pi < +oo is a continuous variable associated to the site i and V{4>i) is an even 
polynomial in the variable In this study, we have only considered the specific model in which 
V{4>i) = (t>i + g{(pi — 1)^, although we can cover interactions of a more general form. 

All these models are expected to belong to the 3D Ising universality class, therefore our eex- 
tensive set of series expansion data can be used to test the accuracy of the basic hypotheses of 
critical scaling and universality with respect to the lattice and the interaction structure, by com- 
paring the estimates of the exponents and of universal combinations of critical amplitudes for the 
various models as well as by forming approximate representations of the equation of state (ES). 
In this study our attitudei^ii^ is, to some extent, complementary to the current one. Usually, uni- 
versality is essentially assumed from the outset: for example, in the renormalization group (RG) 
approach^^"^'', an appropriate scalar field theory in continuum space is taken as the representative 
of the Ising universality class, as suggested by the independence of the renormalization procedure 
from the details of the microscopic interaction. Also in HT and MonteCarlo approaches, attention 
has been recently focusedi^"— on particular continuous- or discrete-spin lattice models which 
exhibit vanishing (or very small) leading non-analytic corrections to scalin g^^i^^ in order to be able 
to estimate more accurately the physical quantities of interest. In this report, we prefer to take 
advantage of our extended expansions to test a wide sample of models, expected to belong to the 
same universality class, and to show how closely, already at the present orders of expansion, each 
model approaches the predicted asymptotic scaling and universality properties. 

The paper is organized as follows: in Sect. II we briefly characterize our expansions, sketch the 
method of derivation and list the numerous tests of correctness passed by the series coefficients. 
In Sect. HI we define the higher-order susceptibilities, whose critical parameters enter into the 
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determination of the scaling equation of state and update an approximate representation of it. In 
Sect. IV we discuss numerical estimates of exponents, amplitudes and universal combinations of 
these, that can be computed from the bivariate series. In the last section, we summarize our results 
and draw some conclusions. 



II. EXTENSIONS OF THE BIVARIATE SERIES EXPANSIONS 

The HT series expansion coefficients for the models under study have been derived by a fully 
computerized algorithm based on the vertex-renormalized linked-cluster (LC) method which calcu- 
lates the mean magnetization per spin in a non-zero magnetic field from the set of all topologically 
distinct, connected, 1-vertex-irreducible (IVI), single-rooted graphsii. We have taken advantage 
of the bipartite structure of the sc and the bcc lattices to restrict the generation of graphs to the 
subset of the bipartite graphs, i.e. to the graphs containing no loops of odd length. 

In the past, the LC method was employed mainly to derive expansions in the absence of mag- 
netic field. In the presence of a field, the most extensive^-^ data so far available in 3D were derived 
indirectly, by transforming^'--^ bivariate LT and high-field expansions'^^ into HT and low-field ex- 
pansions. This computation was performed only for the s = 1/2 Ising model, although some LT 
and high-field data existed also for other values of the spin. Shorter HT expansions in a finite field 
had also been previously obtained^, only for the s = 1/2 model, by a direct expansion of the free 
energy. It is worth noting that we are now in a position to follow the opposite route: namely of 
transforming our bivariate HT data for the spin-s Ising systems into LT and high-field expansions, 
thus extending the known results. 

Our improvements of the presently available HT series in a field are summarized in Table HI in 
the case of the sc and the bcc lattices. Similar extensions for the same class of models, in the case 
of the simple quadratic (sq) lattice and for bipartite lattices in d > 3 space dimensions, will be 
discussed elsewhere. The series expansions coefficients will be tabulated in a separate paper. 

The feasible correctness checks of our computations are inevitably partial, since the extended 
expansions include information much wider than that already available in the literature. The easiest 
non-trivial check is that our procedure yields the known bivariate expansion of the free energy for 
the spin 1/2 Ising model in a finite field on the one-dimensional lattice. We have also checked that 
our results agree, through their common extent, with the old data cited abov^A for the spin 1/2 
Ising system in a magnetic field, both on the sc and the bcc lattices. Otherwise, our results can 
only be compared with the related data in zero field, in particular with the HT expansions of the 
free energy and its second field-derivative, both for the Ising model with general spin-s and for the 
scalar field model, on the sc and the bcc lattices, which have been tabulate d^^i^^ through order 
K'^^, while the 4th field-derivative is already known^^'iS through K^^ for both lattices. Our results 
agree, through their common extent, also with the expansions of the 6th field derivative in zero 
field, tabulated^i up to order K^^, and of the 8th field-derivative, tabulated^^ up to order K^"^ , 
in the case of the spin 1 /2 model on the bcc lattice. We have finally checked that our expansions 
reproduce the sc lattice calculations of the 6th field-derivative (known up to order K^^ ), of the 
8th (known up to order K^'^) and of the 10th (known up to order K^^) in the case of the sc lattice 



scalar field with quartic coupling g = 1.1, which have been tabulated in Ref.jl6|. 

It is fair to remark that the finite-lattice (and the related transfer-matrix) methods of 
expansion2ii2^ have shown more efficient^^ than the LC approach, at least for s = 1/2 in d = 2 di- 
mensions, even in presence of a magnetic field, while they remain rather difficult and unpractical in 
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TABLE I: Maximal order in K of the high-temperature and low-field expansions of the free energy 
for the models in the Ising universality class considered in this note. 





Existing data[4] 


This work 


sc lattice 






Ising s = 1/2 


17 


24 


Ising s > 1/2 





24 







24 


bcc lattice 






Ising 5=1/2 


13 


24 


Ising s > 1/2 





24 







24 



higher space dimensions. In the case of the two-dimensional spin 1/2 Ising model, with the support 
of a variational approximation, these methods made a representation of the ES of unprecedented 
accuracy^'' possible. In the future, these techniques might prove to be competitive^^, in the 3D case 
also for calculations in a nonvanishing field. We believe, however, that it has been worthwhile to 
test and develop also a LC approach, which expresses the series coefficients in terms of polynomials 
in the moments of the single-spin measure and therefore, unlike the finite lattice method, is flexible 
enough to apply also to non-discrete state models such as the one-component scalar field model^^ 
within the Ising universality class studied here and, more generally, to the 0(N)-symmetric spirt^^ 
or lattice-field systems in any space dimension. 

A. The algorithms 

To give an idea of the strong points of our graphical algorithms, we mention that, using only 
an ordinary desktop personal quad-processor-computer with a 4G fast memory (RAM), our code 
can complete in seconds all calculations already documented^*^ in the literature (see Table H]) . The 
whole renormalized calculation presented here, can be completed in a CPU-time of a few days, 
most of which goes in producing the highest order of expansion. In what follows all timings are 
single-core times. 

The LC computation has been split into three parts. First, we generated the simple, bipartite, 
unrooted, topologically distinct IVI graphs. This part is memory intensive, but takes only a few 
hours. Table HI] lists the numbers of these graphs from order 4 through 24. In a second step, we 
computed the single-rooted multigraphs, their symmetry numbers and their lattice embeddings. 
This part of the calculation requires little memory: in the case of the bcc lattice, completing the 
24th order took approximately a day, while in the case of the sc lattice two weeks were necessary. 
In the latter case, most of the time was spent to determine the graph embeddings. These two 
parts of the calculation were implemented by C++ codes, and used the "Nauty"--^ library to 
compute the graph certificates and symmetry factors. The relevant procedures of this package 
were supplemented with the GNU Multiprecision Arithmetics Library^'' to get the exact graph 
symmetry numbers. The third step implements the algebraic vertex-renormalizationii procedure 
by deriving the magnetization from the single-rooted IVI graphs and then, by integration, the free 
energy J^{K,h) = ^ fn{h)K^ . Here K = J/ksT, with ks the Boltzmann constant and T the 
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TABLE II: The number of simple, connected, bipartite, unrooted IVI graphs with / hues and 
with a given number v of odd vertices, which contribute to the HT expansion coefficient of the free 
energy at order 





V 
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10 


12 


14 


16 


totals 
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1 
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10 
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5 




















14 


11 





11 
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2 

















20 


12 
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20 


31 


4 


1 














65 


13 





49 


53 


22 

















124 


14 


20 


101 


194 


54 


7 














376 


15 





258 


432 


238 


20 


2 











950 


16 


84 


520 


1471 


732 


127 














2934 


17 





1482 


3725 


2886 


434 


29 











8556 


18 


300 


3243 


12233 


9531 


2403 


97 


5 








27812 


19 





9646 


33608 


36067 


9675 


845 











89841 


20 


1520 


21859 


109796 


123543 


46241 


4023 


133 








307115 


21 





68697 


318283 


460225 


191416 


26435 


594 


13 





1065663 


22 


8186 


163780 


1048349 


1608030 


858792 


134409 


6672 








3828218 


23 





533569 


3166399 


5970246 


3566324 


757696 


40686 


744 





14035664 


24 


52729 


1328836 


10594514 


21241772 


15475018 


3796365 


317259 


4267 


38 


52810798 



temperature, while h = mH/kBT is the reduced magnetic field. The magnetization is expressed in 
term of the bare vertices M^{h) obtained deriving i times with respect to h the generating function 



Ml^{h) = In 
field system, M^{h) 



sinh(fe(2s+l)/2s) 
sinh(/i/2s) 



in the case of the spin-s Ising systems, or, in the case of the scalar 



In 



-V(<j))+h4> 



For example, the HT expansion coefficient of the free 



energy at order K^, on the hcc lattice, is given by 



while on the sc lattice 



h{h) = 2{Mlf{h)+32{M^)\h)Ml{h) 



h{h) = -iMl)\h) + 18iM^)\h)M^{h). 



(3) 



(4) 



Table IIIII lists the number of monomials of the bare vertices, with a given number v of odd indices 
which contribute to the free energy HT expansion coefficient at order Equivalently, this is the 
number of admissible vertex-degree sequences of the (far more numerous) graphs contributing to 
this coefficient. Notice that the monomials containing at least two bare vertices of odd order are the 
overwhelming majority. They all vanish in zero field, which shows that the finite field calculation 
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TABLE III: The number of monomials in the bare vertices, with a given number v of odd vertices, 
which contribute to the HT expansion coefficient of the free energy on a bipartite lattice, at order 
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o 
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v > 4 
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1 
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11 
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222 
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1379 


2289 


6213 


9881 


18 


278 


1796 


3314 


9566 


14954 


19 





2577 


4635 


14487 


21699 


20 


480 


3370 


6492 


21662 


32004 


21 





4711 


9010 


32134 


45855 


22 


760 


5965 


12430 


46887 


66042 


23 





8257 


16858 


67949 


93064 


24 


1273 


10664 


22895 


97543 


132375 



has a substantially higher complexity. The renormalization through 24th order was performed in 
a few hours. The third step of the calculation is based on codes written in the Python and Sage^ 
languages. 

It is also not without interest that in a preliminary step of our work, we have been able to 
employ the simple unrenormalized linked-cluster method^i, which uses all topologically distinct 
unrooted connected graphs (including multigraphs) to compute the bivariate expansions of the 
free energy through order 20. It takes only one day to complete this calculation. Of course, while 
the unrenormalized procedure is algebraically straightforward, it would make further extensions 
of the series unpractical, using our desktop computers, for the rapid increase with order of the 
combinatorial complexity and, as a consequence, of the memory requirements. The computation 
of the 21st order does not fit in 4 GB of RAM, but would require some increase of memory. These 
calculations are however interesting by themselves, both because the unrenormalized method is 
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still generally (and too pessimistically) dismissed as unwieldy beyond just the first few orders, and 
because they provide a valuable cross check, through order 20, of the results of the algebraically 
more complex vertex-renormalized procedure, which remains necessary to push the calculation to 
higher orders. 

III. ASYMPTOTIC SCALING AND THE EQUATION OF STATE 

The hypothesis of asymptotic scaling^""— for the singular part Fs (t, h) of the reduced specific 
free energy, valid as both h and r approach zero, can be expressed in the form 

Fs{T,h)^\T\^-'^Y±{h/\Tf'). (5) 

where r = (1 — T^/T) is the reduced temperature. The exponent a specifies the divergence of the 
specific heat, (3 describes the small r asymptotic behavior of the spontaneous specific magnetization 
M on the phase boundary (/i — )■ 0"^, r < 0) 

M w B{-Tf (6) 

with B the critical amplitude of M. The exponent 5 characterizes the small h asymptotic behavior 
of the magnetization on the critical isotherm (/i 7^ 0,r = 0), 

\M\^BM^I' (7) 

and Be is the corresponding critical amplitude. For the exponents a and /3, we have assumed 
the values a = 0.110(1) and /3 = 0.3263(4), obtained using the scaling and hyperscaling relations, 
from the HT estimates^^ of the susceptibility exponent 7 = 1.2373(2) and of the correlation-length 
exponent v = 0.6301(2). 

The functions Y±{w) are defined for < w < 00 and have a power-law asymptotic behavior 
as w — )• 00. The + and — subscripts indicate that different functional forms are expected to 
occur for r < and r > 0. The usual scaling laws follow from eq.([5]). The simplest consequence 
of eq.dS]), which will be tested using our HT expansions, is that the critical exponents of the 
successive derivatives of 7s (r, h) with respect to h at zero field, are evenly spaced by the quantity 
A = usually called "gap exponent" . More precisely, let us define the zero-field n-spin connected 
correlation functions at zero wavenumber (also called higher susceptibilities when n > 2) by the 
equation 

Xn{K) = [d^F{h,K)/dh^)h=Q= Yl <SlS2...Sn>c- (8) 

S2,S3,...,S„ 

For odd values of n, these quantities vanish in the symmetric HT phase, while they are nontrivial 
for all n in the broken-symmetry LT phase. For even values of n in the symmetric phase, and for 
all n in the broken phase, scaling implies that, as T — t- along the critical isochore (/i = 0, r > 0) 
or, as T — 7- T~ along the phase boundary, we have 

Xn{r)^Ct\rrH^ + bt\rf + ...) (9) 

where 7n = 7 + ('t- — 2)A, and d are, respectively, the amplitude and the exponent which 
characterize the leading non-analytic correction to asymptotic scaling. The value^'' 6 = 0.52(2) 
has been estimated for the universality class of the 3D Ising model. Assuming also the validity of 
hyperscaling, we can conclude that 2A = 3z^ + 7. 

An important bonus of our bivariate calculations, is the significant extension the HT expansions 
of the higher susceptibilities. We have added one more term to the existingi^ HT expansion of 
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X4:{K), five terms to that-'^ of xei^), seven to thalji^ of XsiK) and nine to tliat of xio{K). In 
the case of the susceptibilities of order 2n > 10, no data at all were available so far. We have 
now extended, uniformly in the order, the HT expansions of all higher susceptibilities X2n{K) with 
2n > 4. In this paper, we shall present only a preliminary analysis of these quantities, while a 
more detailed discussion of our bivariate expansions will be postponed to a forthcoming article. 

The scaling form of the equation of state M = M(h,T), relating the external reduced magnetic 
field h, the reduced temperature r and the magnetization M, when h and r approach zero, is 
simply obtained by differentiating eq.([5]) for /s(r, h) with respect to h 

M ^ -\T\I^Yi^\h/\T\l^^) (10) 

Here we have used the relation 7 = — 1). By further differentiation of eq. (jlOp with respect to 
the field, also the higher susceptibilities are recognized to have a scaling form 

Xn{h,T) = {d^-^M/dh^-^) « -\T\-^-Yt\h/\T\P'). (11) 

The hypothesis of universality states that, in addition to the critical exponents, the function Y±{w), 
(and therefore also its n—th derivative Y^\w)) is universal^ up to multiplicative constants {metric 
factors^^) which fix the scales of h and r in each particular model within a universality class. 
Accordingly, one can conclude that a variety of dimensionless combinations of critical amplitudes 
are universal. 

The ES can also be written in the equivalent for m^^'^^'^*^ 

h{M, r) w M\M\^'^f{T/\M\^/^) (12) 

in which a single scaling function f{x), universal up to metric factors, describes both the regions 
r < and r > 0. The function h{M, r) is known^- to be regular analytic in a neighborhood of 
the critical isotherm and of the critical isochore. From general thermodynamic arguments.^2. one 
can infer that f{x) is a positive monotonically increasing regular function of its argument, in some 
interval —xq < x < 00, with xq > 0. Moreover /(— xq) = 0. The local behavior of the function 
f{x) can be further determined, by the requirement of consistency with the scaling laws, in terms 
of critical amplitudes of quantities computable from our HT and LT series. By differentiating this 
form of the ES with respect to M, we get the asymptotic behavior /(x) oc x'^ , for large positive x. 
Setting r = 0, the ES reduces to eq.(I7|) and /(O) = . If /i — )• at fixed r < 0, we expect to 
find a nonvanishing spontaneous magnetization M, therefore the ES implies that /(x) must vanish. 
Since /(— xq) = 0, we have —t/M^/^ = xq and, from eq.®, we conclude that xq = B^^^f^. We 
can then fix the metric factors by normalizing the field to B^^ and the reduced temperature to 
The expansion of /(x) for large positive x is expressed in terms of the critical parameters 
characterizing the HT side of the critical point. The small x expansion, which uses the parameters 
of the critical isotherm, and the negative x region related to the parameters of the LT side of Tc, 
will be discussed in a forthcoming paper presenting our analysis of the extended LT expansions. 



Summarizing the more detailed discussion of Ref.[2l|], we can also observe that, in the large pos- 
itive X (small magnetization) region, where the magnetic field h{M, r) has a convergent expansion 
in odd powers of M, the ES is more conveniently expressed in terms of the variable z = Mt~^Xq. 
The ES takes then the form 

h{M,T) = h\Tf^F{z) (13) 

where his a constant and F{z) is normalized by the equation -F'(O) = 1. The small z expansion of 
Fiz) can be written as 
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F{z) = z + -z^ + F5Z^ + F7z'^ + ... (14) 

The coefficients F^,F-j,... are defined by the equation F2n-i = f'tn/i'^''^ ~ ™ terms of 
the ratios which will be introduced in the next section. They have been computed within 
the RG approach^^i^, by the e-expansion (e = 4 — d) up to five loops, by the perturbative g- 
expansion at fixed dimension d = 3 up to the same order, by other RG approximations^^^-^I, by 
HT expansion o^^i^^ , by MonteCarlo methods^*^. Our estimates of the first few by extended 
HT expansions, are reported in Table HXl 



A. A parametric form of the ES 

A parametric fornt^""— has been introduced to formulate an approximate representation of the 
ES in the whole critical region and as an aid in the comparison with the experimental data. The 
parametrization is chosen to embody the analyticity properties of /i(M, r) and the scaling laws. 
These properties make the parametric form convenient to approximate the ES in the whole critical 
region by using only an HT input, such as the small z expansion eq. (|14p of F{z). In this approach, 
the scaled field and the reduced temperature are expressed as the following functions 

M = moR^e (15) 
T = R{1- 9^) (16) 

h = hoR'^^l{e) (17) 

of generalized radial and angular coordinates R> and —9q<9<9q, with > 1 the smallest 
positive zero of the function l{9). The radial coordinate R measures the distance in the /i, T plane 
from the critical point, and the angular coordinate 9 specifies a direction in this plane. Therefore, 
= corresponds to the critical isochore, = ±1 is associated to the critical isotherm and 9 = ±^o 
to the coexistence curve. The function l{9), normalized by /'(O) = 1, is odd and regular for |0| < ^o, 
as implied by the regularity of /(x) and the invertibility of the above variable transformation in 
this interval. 

The variable z is then expressed as 

(18) 



{l-9^Y 

and the function F{z) of eq. lfH]) is related to l{9) by 

l{9) = -^{{l-9y+^F{z{9)) (19) 

Here p = mo^Q is a positive constant related to the arbitrary normalization constant mo ap- 
pearing in eq. (llSh . If F{z) were exactly known, the corresponding l{9) given by eq. (I19p should 
not depend on p. However the polynomial truncations of l{9) that can be formed from the first 
few available terms of the expansion eq. of F{z) will have coefficients l2n+i{p) depending not 
only on the coefficients F^,Fi, ... and on the exponents f3 and 7, but also on p. 
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In particular—, expanding both sides of eq.([T9]), one obtains: 

h{p) = \p'-l (20) 
kip) = ^7(7 - 1) + ^(2/3 - 7)/,2 + F5/ (21) 

Hp) = - 1)(7 - 2) + ^(2/3 - 7)(2/3 - 7 + l)/,^ + (4^ _ ^)F5/ + Fjp^ (22) 

«9(/3) = ^7(7 - 1)(7 - 2)(7 - 3) + ^(2/3 - 7)(2/3 - 7 + l)(2/3 - 7 + 2)/>2 + (23) 

^(4/3 - 7) (4/? - 7 + 1)^5/ + (6/3 - 7)F7p6 + ^g^s 

etc. 

The dependence on p of the coefficients ^2™+! has been exploited to improve the approximation 
of l{9). A first approach consists in fixing p to the value pm which minimizes^ the modulus of 
the highest-order expansion coefficient l2n+i{p) of l{9) that can be determined reliably from the 
available coefficients i<2n-i- A second methodic is based on computing some universal combinations 
of critical amplitudes in terms of l{6) and then in choosing for p the unique value that makes all 
such quantities stationary. We may follow this route and consider, for example, the dependence 
on p of the universal ratio of the susceptibility amplitudes above and beneath Tc, namely C^/C^, 
and of the ratios C^B"^/ (€2)^ and C2 B^^^ / B^. If we plot these quantities vs p^, we obtain FiglU 
which indicates that the choice p^ = 2.615 should be optimal. Then, using the central values both 
of the coefficients F2n-i up to n = 7, as obtained from our Table HXl and of the exponents (5 and 7 
as indicated above and fixing p to its optimal value, the following form of l{9) can be determined 

l{e) f^e- 0.8014(50)6l3 + 0.00946(30)6*^ + 0.00141(40)6*^ + 0.00029(10)6*^ - 0.00011(5)6*^^ (24) 

Here we have neglected the term in 6^^, whose coefficient is 0(10^^), and have indicated the last 
three terms only to show that their contribution in the interval of interest \6\ < Oq is very small. 
The function 1(9) vanishes at 6 = Oq ^ ±1.1273. 



The analogous result for this auxiliary function obtained in Ref . |22l| . fixing p by the first method 



and choosing the values /3 = 0.3258(14) and 7 = 1.2396(13) of the critical exponents, is 

1(6) ^e- 0.762(3)0^ + 0.0082(10)0^ (25) 

which vanishes at ~ ±1.1537. In this case, the coefficients -F2n-i were obtained by a RG five-loop 
perturbation expansion in d = 3. On the other hand, computing the -F2n-i by the RG e-expansion 
to fifth order and choosing /3 = 0.3257(25) and 7 = 1.2355(50) leads^^ to: 

1(6) ^e- 0.72(6)0^ + 0.0136(20)0^ (26) 

More recently in Ref. [16], using values of the exponents very near to those used in our paper, 
and deriving the -F2n-i from an HT expansion of sc-lattice scalar-field models with self-couplings 
appropriately chosen to suppress the leading correction to scaling, the following expression was 
obtained, 

1(6) ^e- 0.7367430^ ± 0.0089040^ - 0.0004720^ (27) 



which vanishes at ~ ±1.1741. As stressed in Refs. 2l|j22l ]. the alternative forms eqs. (I24p . (I25p . 
(j26|) and (p7|) cannot be directly compared, because they are associated to different parametriza- 
tions (different values of p). One should rather compare the universal predictions obtained from 
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TABLE IV: A few universal amplitude combinations obtained in this work from the parametric 
form eq. (|24p of the ES. For comparison, we have reported also the results obtained from: the 
slightly different parametric form eq. (j27D of the ES in Ref. la] based on shorter HT expansions 
and from the parametric forms in Ref. [2^ obtained, either from the RG e— expansion or from 
the 5— expansion and using values of the exponents slightly different from those used to get the 
estimates in the first two columns. 



Universal ratios 


This work 


Ref. [16] 


e— expans.[22] 


5— expans.[22] 




0.530(3) 


0.529(6) 


0.527(37) 


0.537(19) 




4.78(3) 


4.78(5) 


4.73(16) 


4.79(10) 




-9.2(3) 


-9.3(5) 


-8.6(1.5) 


-9.1(6) 


-CtB^/{Ct? 


7.8(1) 


7.83(4) 


8.24(34) 




C+B'-^/Bi 


1.66(2) 


1.665(10) 


1.648(36) 


1.669(18) 


aA+Ct/B^ 


0.0563(5) 


0.0562(1) 


0.0569(35) 


0.0574(20) 


-C^B/iC,)^ 


6.015(15) 


6.018(20) 


6.07(19) 


6.08(6) 



them, for example, for the universal amplitude combinations, which we have reported in Table 
IIVI In this Table, A~^ and A~ denote the amplitudes of the specific heat above and beneath Tc. 
Our estimates are compared with the corresponding ones obtained^^ from the polynomials l{9) in 
eq. (|25p and eq. (j26p and with those obtainedi^ from eq. (j27p . The results from the various methods 
show a good overall consistency. In Table IIVI we have not reported MonteCarlo estimates, which 
also are available^i^i^ for the ratios A~^/A~, C^i?^~^/S^ and C^/C^. It is worth to remark 
that the computation of the first quantity is difficult, because the weak singularity of the specific 
heat forces to extend the simulation very close to the critical point. Recent estimates^*^ of this 
ratio, in the range 0.532(7) — 0.540(4), based on simulations of large lattices, might now super- 
sede older results, which were ~ 6% larger, thus improving the agreement with the ES estimates 
of Table ITVl Also the second ratio, involving the amplitude Be of eq.© is difficult to compute 
by simulations, for similar reasons. The result C^^B^-^/B^ = 1.723(13) of Ref. [29], is somewhat 
larger than the estimates from the ES. In the case of the third ratio, recent simulations^^'^^ have 
changed the previous larger estimates to values in the range 4.67(3) — 4.72(11), closer to the ES 
results. In the Table IIVI we also have not shown the few available experimental estimates of some 
of these combinations, nor those, tabulated in Refs. 16lj2l|j22l ]. which are based on direct evalua- 



tions of the amplitudes by LT and HT expansions. They are completely compatible with the ES 
results of Table HVj but the comparison is not stringent, due to the large uncertainties. Therefore 
we plan to improve the series determinations of the amplitudes on the critical isotherm and on 
the coexistence curve, exploiting our extended bivariate LT expansions of the free energy for the 
spin-s Ising models. A more detailed analysis of our results, along with estimates of other universal 
amplitude combinations, and a wider comparison among the results in the literature is deferred to 
a forthcoming paper. 



B. The HT zero-momentum renormalized couplings 

The HT expansions of the higher susceptibilities^ will be used to evaluate their critical ampli- 
tudes and correspondingly the critical limits of the zero-momentum n— spin dimensionless HT 
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renormalized couplings (RCC's). These quantities enter into the approximate forms of the scahng 
ES eqs. (fTO]l and (fT2|). 

In the HT phase the first few 2n-spins RCC's are defined as the critical limits as K ^ of 
the following expressions 

V X4iK) 



9t{K) 



9tm 



9t{K) 

y2 



X8{K) ^^^X6{K)X4{K) 



xUk) 



280(:^)^ 

^xUk)J 



^,,,X.{K)xm ^,,,XI{K) 



xliK) 



xliK) 



-4620 



{K)) 



(28) 
(29) 
(30) 

(31) 



y5 



Xi2{K) _^^^qXio{K)x4{K) 



XliK) 



xliK) 



792 



X8iK)x6{K) 
xliK) 



.,j,,0 XsiK)xliK) _ 3gQ3gXlWx4W ^ ,Qo,,o XeiK)xliK) _ ^.^..^^(X^y 
XUK) xUK) xliK) ^ xliK) J 



(32) 



9UK) 



e^K) L xliK) 



X.iK) ^3^^Xi.(i^)X4(i^)_3^^,^Xio(5)xl(K) ^33^ 



xliK) 



xliK) 



+2002^i°i^M^ + 17164^^ + 3203200^^^^^MW _ ^Q^^Q^XsiK)xeiK)xm 



xliK) 



-126126 



xliK) 



+ 10090080 



xliK) 
xUK)xliK) 



xfiK) 



xliK) 



95295200^^^^^?M^ + 190590400 f^^^V 
X2iK) ^xiiK)J 



■xl{K) xfiK) 
Here ^iK) is the second moment correlation-length, defined by 

6X2 

with ^2 the second moment of the correlation function expressed as 

/X2 = < Sofia; >c . 



(34) 



(35) 



Both the HT expansions of xiK) and ^2iK) are tabulatedi^ through order K^^ for the spin-s Ising 
models. 

The volume V per lattice site takes the value 1 for the sc lattice and 4/3^/3 for the bcc lattice. 
The definitions of the quantities (7^(-fC) given here differ by a factor (2n)! from those of Ref. 51 1. 
Also the quantities, 

X^{K)x2n+AiK) 



I 



2n+4 



(K) 



XT\K) 



(36) 



with n > 1, whose critical values are the universal amplitude combinations first described^® in the 
literature, and the closely related quantities 

92niK) 



r2niK) 



94iK) 



n-l 



(37) 
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which share the computational advantage of being independent of the correlation length, will be 
of relevance in what follows. The finite critical limits g^, and /^^4 of the RCC's, of the ratios 
r2n{K) and of the quantities l2n+i{K), represent universal combinations of HT amplitudes that 
should be considered together with those listed in Table llVi We have not included the expressions 
of higher-order RCC's, because, in spite of our extensions, the available series might not yet be long 
enough to determine safely their critical limits. One should notice that, from the point of view of 
numerical approximation, the and also the quantities derived from them like r^, are difficult 
to compute, unless 2n is small, because they result from relatively small differences between large 
numbers. These estimates can be reliable provided that the uncertainties of the large numbers 
are much smaller than their difference. For the same reason, these quantities are notoriously even 
more difficult to compute by stochastic methods. 

IV. METHODS AND RESULTS OF THE SERIES ANALYSIS 
A. Extrapolation methods 

In the numerical analysis of the series expansions of physical quantities, we shall follow two 
procedures aimed to determine the critical parameters, namely the values of these quantities at 
the critical point, whenever they are finite, or if they are singular, the locations, amplitudes and 
exponents of the critical singularities on (or nearby) the convergence disk in the complex K plane. 

A first procedure used in our series analysis, is the differential approximant (DA) method^, a 
generalization of the well known Fade approximant method^^, having a wider range of application. 
In this approach, the values of the quantities or the parameters of the singularities can be estimated 
from the solution, called differential approximant, of an initial value problem for an appropriate 
ordinary linear (first- or higher-order) inhomogeneous differential equation. This equation has 
polynomial coefficients defined in such a way that the series expansion coefficients of its solution 
equal, up to a certain order, those of the series under study. The various possible equations, and 
therefore the various DAs that can be formed by this prescription, are usually identified by the 
sequence of the degrees of the polynomial coefficients of the equation. The approximants are called 
first-order, second-order DAs etc., according to the order of the defining equation. The convergence 
of the procedure in the case of the Ising models can be improved by first performing in the series 
expansions the variable transformation^^ 



aimed at reducing the infiuence of the leading corrections to scaling. Here 6 is the exponent 
which characterizes these corrections. A sample of estimates of the parameters of the critical 
singularity is obtained from the computation of many "quasi diagonal" DAs, namely approximants 
with small differences among the degrees of the polynomial coefficients of the defining differential 
equation, which use all or most of the given series coefficients. A first estimate of a parameter, 
along with its uncertainty, results from computing the sample average and standard deviation. The 
result can then be improved by discarding from the sample single estimates which appear to be 
obvious outliers, and recomputing the average of the reduced sample. A conventional guess of the 
uncertainty of the parameter estimate is finally obtained simply, and rather roughly, as a small 
multiple of the spread of the reduced sample around its mean value. This subjective prescription 
might, to some extent, allow for the difficulty to infer possible systematic errors, and to extrapolate 
reliably a possible residual dependence of the estimate on the maximum order of the available series. 




(38) 
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A second approach is based on a faster converging modification of the standard analysis of 
ratio-sequence of the series coefficients and will be denoted here as the modified ratio approximant 
(MRA)^^ technique. Let us assume that the singularity of the series expansion of a physical 
quantity, which is nearest to the origin of the complex K plane, is the critical singularity, located 
at Kc and characterized by the critical exponent A and the exponent 9 of the leading correction 
to scaling, (this hypothesis is generally not satisfied for the LT series). Then eq.([9]) implies the 
following large r behavior of the series coefficients cv 

"^-1 '- r(A) b 



C 



r(A) 



k: 



1 + 



+ 0(l/r) 



(39) 



r(A - 9) rO 

In this case, the MRA method evaluates Kc by estimating the large r limit of the approximant 
sequence 

' Cr-2Cr-.3V/'^ [ Sj- + Sr-2 



Cf' — 1 



exp 



with 



In 



-1 



+ In 



-1 



/2. 



(40) 



(41) 



By using the asymptotic form Eq. (j39p . we can obtain^- the large r asymptotic behavior of the 
sequence of MRA approximants of the critical inverse temperature 

r(A) 9'^{l - 



Kr 



1 



+ 0(l/r" 



2r(A - 9) ri+' 

The method estimates also the critical exponent A from the sequence 

{Sr + Sr-2) 



(42) 



(A)r 



1 + 2- 



{Sr - Sr-2 

with Sr defined by eq. ([^T]) . In this case the large r asymptotic behavior of the sequence (A)^ is 

r(A) 9{l-9'^)b 



(43) 



(A). 



A 



r(A 



+ 0(l/r) 



(44) 



If the available series expansions are sufficiently long (how long cannot unfortunately be decided o 
priori), the estimates of the critical points and exponents obtained from extrapolations based on 
eqs. (j42p and (|44p can be competitive in precision with those from DAs. If, on the other hand, 
the series are only moderately long or the exponent A >> 1, then corrections of order higher than 
1/r^^^ in eq. (|42p (or higher than 1/r^ in eq. (|44p ) might still be non-negligible. The same remark 



applies if the 0(l/r) terms in eq. p9p are not sufficiently small. Therefore, in some cases, eqs. (|42p 
and (|44p might be inadequate to extrapolate the behavior of the few highest-order terms of the 
MRA sequences. 



B. Critical parameters of the higher susceptibilities 

For both methods sketched in the previous paragraph, the main difficulties of the numerical 
analysis of the HT expansions are related to the presence of the leading non-analytic corrections 
to scaling which appear in the near-critical asymptotic forms of all physical quantities. It was 
however observed^ft!^- that the amplitudes of these corrections are non-universal and therefore, by 
studying families of models expected to belong to the same universality class, one might be able 
to single out special models for which these amplitudes have a very small or vanishing size. These 



14 



models would then be good candidates for a high-accuracy determination of the critical parameters 
of interest. In the literature, various models which share this property to a good approximation, 
have been subjected to analysis: among them, the lattice cp'^ model on the sc lattice with the value 
g = 1.1 of the quartic self-couplin g^^ '^^ , or the same model on the bcc lattice with the coupling^- 
g = 1.85. Also the spin-s = 1 and s = 3/2 Ising systems on the bcc lattice^^, show very small 
corrections to scaling. All these models will be considered here. 

An accurate estimate 2A = 3.1276(8) of the gap exponent which improves the four-decade-old^ 
estimate 2 A = 3.126(6), based on 12th order series, had been already obtained from the known 
23rd order HT expansions of Xi{K) for the spin-s Ising modelsiS and for the lattice scalar field^-, 
on the sc and bcc lattices. The addition of a single coefficient to the expansion of x^i^) does not 
urge resuming a full discussion of the estimates of this exponent and of the validity of hyperscaling 
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on the HT side of the critical point, already tested with good precision in Refs. 

To get some feeling of the reliability of the estimates that can be obtained from a study of our 
HT expansions of the higher susceptibilities X2niK), it is convenient to test how accurately the 
critical inverse temperature Kc and the critical exponents 72n can be determined from them by 
using MRAs. Let us for example consider the above mentioned self-interacting lattice scalar field 
model of eq.([2]) on the sc lattice with quartic coupling g = 1.1. In Fig. [21 we have plotted vs r^^^ 
the sequences of the MRA estimates {Kc)r for Kc, as obtained from the HT expansions of X2n{K) 
with 2n = 2,4, ...22. The MRA sequences are normalized by the appropriate limiting value of the 
sequence {Kc)r, estimated in Ref.[l8|] and reported in Table IVl to make them easily comparable 
with the corresponding sequences obtained from other models in the same universality class. The 
choice of the plotting variable is suggested by eq. (|12|) . For the susceptibilities of order 2n > 6 the 
curves indicate the presence of strong corrections 0(r~'^), with a between 3 and 5, and show that 
simply using eq. (jl2|) . at the present orders of expansion, would be inadequate for extrapolating to 
r — )• oo the MRA sequences. No significant quantitative difference in behavior is observed in the 
analogous plots for the other models examined in this study, even for those with non-negligible 
amplitudes of the leading corrections to scaling. On the contrary, in other cases, for example for 
the Ising model with spin-s = l/2ors = lon the bcc lattice, the convergence looks even slightly 
faster. From these plots one may conclude that, as the order 2nof the susceptibility X2n{K) grows, 
increasingly long expansions are needed^ in order that the MRA sequences reach the asymptotic 
form eq. (|42|) and therefore a given precision can be achieved in the estimate of Kc- The general 
features of this behavior can be tentatively explained arguing""- that the dominant contributions 
to the HT expansion of X2niK), at a given order K^ , come from those spin correlation functions 
in the sum of eq.®, for which the average distance among the spins is ~ r/2n. Accordingly, it 
seems that the presently available expansions of the quantities X2n (K) , in spite of having the same 
number of coefficients, might not have the same "effective length" , because they describe systems 
which are, in some sense, rather "small", the more so the larger is 2n . One might conclude that 
the estimates of the critical quantities, derived from the X2n{K), should probably be taken with 
some caution for large n, even in the case of models with very small leading corrections to scaling. 
However, in what follows, we shall observe that, in some cases, in spite of these difficulties, the 
DAs seem to yield smooth and reasonable extrapolations of these series to the critical point. 
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TABLE V: Estimates [iTlllSl ] of the critical inverse-temperatures Kc used in our study of the Ising 
systems with spin s and of the lattice scalar field systems, on the sc and the bcc lattices. 





s = 1/2 


s = 1 


s = 3/2 


s = 2 


s = 5/2 


s = 3 




jy-bcc 


0.221655(2) 
0.1573725(10) 


0.312867(2) 
0224656(1) 


0.368657(2) 
0.265641(1) 


0.406352(3) 
0.293255(2) 


0.433532(3) 
0.313130(2) 


0.454060(3) 
0.328119(2) 


0.375097(1) 
0.2441357(5) 



C. Scaling and the gap exponent 

For Ising models on the bcc lattice with spin s = 1/2, 1, 3, we have computed the sequences 
of estimates of the exponent differences Z)„ = j2n — l2n-2, with n = 2,3, ..11. These estimates 
are obtained from second-order DAs of the ratios X2n{K) / X2n-2{K) , which use at least 19 series 
coefficients. We have imposed that the critical inverse temperatures, for the various spin systems, 
take the appropriate valuesi^, listed in Table El In Figl3l the exponent differences D„ are plotted 
vs n. We have observed above that, as a simple consequence of the scaling hypothesis, when the 
maximum order of the available HT series grows large, the Dn should all converge to the same 
value, equal to twice the gap exponent A, thus being independent of the order 2n of the higher 
susceptibilities entering into the calculation. For some particular values of the spin, e.g. s = 1 and 
s = 3/2, our central estimates depart by less than 0.1% from the expected result, for all values of 
n considered here. For other values of the spin, e.g. s = 1/2, a residual spread of the data remains, 
which however is quite compatible with the errors due to the finite length of the series and to the 
likely presence of sizable corrections to scaling, particularly when n is large. This computation 
can be repeated, with similar results, but somewhat larger error bars, for the spin-s Ising system 
on the sc lattice. We have shown in Figj4] the results of the same computation for the two lattice 
scalar field systems with suppressed leading corrections to scaling, studied here. In view of the 
above remarks concerning the effective length of the expansions of the X2n [K) , our results confirm 
the expectation that, in general, the uncertainties of the results should grow with n. It should be 
noted that, both in the case of the spin-s Ising system and of the scalar field system, the bcc lattice 
expansions have a distinctly smoother and more convergent behavior than for the sc lattice, on 
a wider range of values of the order 2n of the susceptibilities, probably because the coordination 
number of the hcc lattice is larger. In conclusion, our results support the validity of the scaling 
property, while the rather accurate independence of the estimates of the gap exponent on the lattice 
structure and, in the case of the Ising models, on the value s of the spin, is a valuable indication 
of universality. Finally, it is worth to stress that the results of Fig. [3] and Fig. d] for n > 4 would 
be difficult to obtain by numerical approaches other than series expansions. 

D. The ratios and the critical amplitudes of the higher susceptibilities 

In Ref.Qj the critical amplitudes of the higher susceptibilities were estimated from 12th 
order series, for the simple s = 1/2 Ising model, assuming the now outdated values 7 = 5/4 
and A = 25/16 for the exponents. These series were not long enough that any estimate of the 
uncertainties could be tried. It is then worthwhile to update the estimates of these amplitudes by 
using our longer expansions and biasing the extrapolations by the more precise modern estimates of 
the exponents and the critical temperatures cited above. We can moreover obtain the corresponding 
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TABLE VI: Our final estimates, by first-order DAs, of the critical amplitudes, of the second- 
moment correlation length eq. ([3l]l , and of the susceptibilities X2n (K) eq. ([9]) , on the HT side 
of the critical point, for Ising models with various values s of the spin on the sc and the bcc lattices 
and for the lattice scalar field with cj)^ self-interaction. The quartic coupling has the value g = 1.1 
for the sc lattice, while g = 1.85 for the bcc lattice. For convenience, we have reported the value 
of C2„/(2n)!. 





^ I/O 

-S — 1 / z 


^ 1 

-S — 1 


s — o/z 


s — z 


5 — 0/z 


5 — O 


w!>4 
H> 


bcc 
















f+ 

h 


U.4Doi(o ) 




C\ A'X C\'7 ( 0\ 

U.4iU/(z ) 


U.4(J4o(i ) 


U.4UlU(i j 




n A'X ACi.( '\\ 

U.4i4o(i j 


/2! 


0.5202(9) 


0.3105(6) 


0.2481(5) 


0.2186(5) 


0.2019(5) 


0.1910(4) 


0.2741(7) 


C7+/4! 


-0.1416(9) 


-0.0377(1) 


-0.02175(7) 


-0.0161(1) 


-0.0134(1) 


f\ f\-\ -t cif\fci\ 

-0.01180(8) 


-0.02728(8) 




0.1224(9) 


r\ r\'i A f f / o\ 

0.01455(8) 


r\ r\r\r'f\r' / A\ 

0.00605(4) 


0.00377(6) 


0.00282(4) 


0.00231(3) 


0.00862(8) 


Cg+/8! 


-0.150(3) 


-0.00798(9) 


-0.00240(2) 


-0.00125(3) 


-0.000845(9) 


-0.000646(7) 


-0.00387(3) 


/^+ /I nl 

<-io/10! 


0.22(i) 


0.0052(1) 


(\ f\f\-\ 1 O / 1 \ 

0.00113(1) 


f\ r\r\c\ A Ci'T ( f\\ 

0.000497(9) 


0.000301(5) 


0.000215(3) 


0.00207(4) 


C+/12! 


-0.35(5) 


-0.0037(1) 


-0.00059(2) 


-0.00022(2) 


-0.000119(5) 


-0.000079(3) 


-0.00123(4) 


C+/14! 


0.57(9) 


0.0029(2) 


0.00032(3) 


0.000101(8) 


0.000050(3) 


0.000031(2) 


0.00077(5) 


sc 
















n 


0.5070(5) 


0.4588(4) 


0.4429(4) 


0.4356(4) 


0.4317(5) 


0.4294(5) 


0.4151(1) 


Ct/2\ 


0.5608(9) 


0.338(2) 


0.270(2) 


0.239(1) 


0.220(1) 


0.208(1) 


0.2384(7) 


C+/4! 


-0.1608(5) 


-0.0432(2) 


-0.0249(2) 


-0.01847(9) 


-0.0153(1) 


-0.0135(1) 


-0.01595(3) 


C+/6! 


0.146(3) 


0.0175(2) 


0.00729(4) 


0.00454(3) 


0.00339(1) 


0.00277(2) 


0.00339(1) 


C+/8! 


-0.187(9) 


-0.0101(2) 


-0.00302(5) 


-0.00158(3) 


-0.00106(2) 


-0.000809(9) 


-0.00102(1) 


ctjm 


0.26(6) 


0.0069(3) 


0.00148(6) 


0.000655(9) 


0.000393(9) 


0.000279(9) 


0.000367(6) 


C+/12! 


-0.19(9) 


-0.0049(9) 


-0.00079(9) 


-0.00030(3) 


-0.000162(8) 


-0.000107(6) 


-0.000146(5) 


C+/14! 


0.026(9) 


0.0016(9) 


0.00034(9) 


0.00013(4) 


0.000066(9) 


0.000041(8) 


0.000061(4) 



TABLE VII: Estimates of the amplitudes eq. ([9]), tabulated in Ref.fa] without indication of 
error and only in the case of the Ising model with s = 1/2. 





C+/2\ 


C+/4! 


C^/6\ 


Ct/8l 


C+/10! 


C+/12! 


sc lattice 


0.5299 


-0.1530 


0.1366 


-0.1722 


0.2601 


-0.4526 


bcc lattice 


0.4952(5) 


-0.1385 


0.1169 


-0.1397 


0.2023 


-0.3297 



informations also for the other models under scrutiny. The critical amplitudes CJ^ with n > 2 
can also be evaluated, with results consistent within their errors , in terms of C^, and of the 
universal critical values 7^_|_4 of the quantities defined by eq. ([36]) . In this approach only the 
estimates of and need to be biased with both the critical temperatures and the exponents, 
while of course the estimates of l2n+4 ^^'^^ to be biased only with the critical temperatures. Our 
final estimates for the amplitudes are collected in the Table IVIi The results of Ref . |^ , are 
reproduced for comparison in Table IVIII 

We have estimated the critical values of the HT expansions of the RCC's either directly, by 
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TABLE VIII: Our final estimates of the universal critical values of the quantities r2n{K), 
with 3 < n < 7, for Ising models with various values s of the spin and for the scalar field with 
(/)^ self- interaction, on the sc and the hcc lattices. The quartic coupling has the value g = 1.1 for 
the sc lattice, while g = 1.85 for the hcc lattice. We have also reported the values of the universal 
amplitude ratios /^_,_4/(2n + 2)!, with 1 < n < 5. 



hcc lattice 


s = 1/2 


s = 1 


s = 3/2 


s = 2 


s = 5/2 


s = 3 




at 


23.56(4) 


23.54(3) 


23.54(3) 


23.54(3) 


23.54(2) 


23.54(2) 


23.56(1) 


r+ 


2.064(8) 


2.061(9) 


2.063(5) 


2.062(5) 


2.062(5) 


2.062(5) 


2.061(2) 


1 


2.54(5) 


2.64(5) 


2.61(4) 


2.60(4) 


2.59(4) 


2.58(5) 


2.54(4) 


1 


-15.1(9) 


-15.0(5) 


-15.4(6) 


-15.9(7) 


-16.1(8) 


-16.3(8) 


-15.2(2) 




45(7) 


40(5). 


44(5) 


48(6) 


51(7) 


53(8) 


44(3) 


1 


1504(240) 


1490(115) 


1359(82) 


1366(100) 


1319(100) 


1300(100) 


1615(120) 


l6/4! 


0.3307(3) 


0.3308(5) 


0.3307(2) 


0.3307(2) 


0.3307(2) 


0.3308(2) 


0.3308(1) 


Ig/e! 


0.2319(6) 


0.2319(5) 


0.2320(4) 


0.2320(4) 


0.2320(4) 


0.2320(4) 


0.2321(2) 


/io/8! 


0.1667(2) 


0.1668(3) 


0.1666(2) 


0.1664(3) 


0.1666(4) 


0.1669(5) 


0.1667(2) 


/12/IO! 


0.1216(3) 


0.1215(1) 


0.1216(2) 


0.1214(2) 


0.1214(3) 


0.1215(2) 


0.1213(3) 


I14/I2! 


0.0894(2) 


0.0892(2) 


0.0894(3) 


0.0894(3) 


0.0893(3) 


0.0893(3) 


0.0897(6) 


sc lattice 
















at 


23.59(4) 


23.57(6) 


23.56(2) 


23.55(4) 


23.55(4) 


23.55(3) 


23.55(3) 


r+ 


2.067(11) 


2.066(8) 


2.064(7) 


2.065(7) 


2.066(7) 


2.065(7) 


2.057(3) 


r+ 


2.51(7) 


2.41(5) 


2.45(10) 


2.57(10) 


2.61(9) 


2.61(9) 


2.45(5) 




-17(2) 


-14(2) 


-14(2) 


-14(1) 


-14(1) 


-14(1) 


-15.4(2) 




45(8) 


44(8) 


44(6) 


52(4) 


54(5) 


51(6) 


62(3) 


'14 


1460(240) 


1390(130) 


1644(115) 


1477(120) 


1362(150) 


1310(150) 


1176(140) 


/6/4! 


0.3306(5) 


0.3306(3) 


0.3307(3) 


0.3306(3) 


0.3306(3) 


0.3306(3) 


0.3310(1) 


4/6! 


0.2320(12) 


0.2316(7) 


0.2316(8) 


0.2317(8) 


0.2318(8) 


0.2318(7) 


0.2324(3) 


/10/8! 


0.1678(8) 


0.1670(10) 


0.1665(4) 


0.1665(3) 


0.1665(3) 


0.1665(3) 


0.1665(5) 


/12/IO! 


0.1211(6) 


0.1214(5) 


0.1213(5) 


0.1210(5) 


0.1209(6) 


0.1208(4) 


0.1206(9) 


/14/I2! 


0.0908(12) 


0.0906(12) 


0.0894(10) 


0.0892(10) 


0.0889(12) 


0.0893(8) 


0.0884(7) 



extrapolatiouiaii^ to of the simple auxiliary function 

W2n{K) = {K/K,f-^g+^{K) (45) 

designed to be regular at K = and therefore more convenient to study by DA's, or, more 
conveniently, but with consistent results, from the computation of the quantities using eq. 
(j37p . In Fig. [SI our estimates for gt are plotted vs the value s of the spin for Ising systems on the 
sc and the hcc lattices and compared to the estimate gt = 23.56(3) (dashed line) of Ref.[l3|. In 
the same figure, we have also shown the values of gt for the scalar field model on both lattices. 

Table IVIIII lists our estimates of the quantities ltn+4 ^^id , obtained from first- and second- 
order DA's, for a few spin-s Ising systems and for the lattice scalar field, on the sc and the 
hcc lattices. We have imposed that the critical inverse temperatures take the appropriate values 
reported in Table |V] and that an antiferromagnetic singularity is present at —Kc- Only for the 
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TABLE IX: Our final estimates of tlie quantities gj, rg , rg , rj^, rj^ and rj^, obtained either from 
the 4>'^ results on the bcc lattice or from a weighted average of the results on both the sc and the 
bcc lattices, are compared to estimates in the recent literature. These have been obtained: i) from 
HT expansions'^, shorter than those analyzed here, of the sc lattice scalar field with (p'^ or (p^ self- 
interactions and appropriate couplings; ii) from the expansion^^ in powers of e = d — 4, within the 
RG approach; iii) from the g— expansio n^-?'"" in fixed dimension d = 3, within the renormalization- 
group approach; iv) from various approximations^^"— of the renormalization-group equations; v) 
from MonteCarlo simulations^!^. 







-6^ 


r+ 




rt 


rt. 


This work 


23.56(1) 


2.061(2) 


2.54(4) 


-15.2(4) 


45(5) 


1400(200) 


HT scalar sc[16] 


23.56(2) 


2.056(5) 


2.3(1) 


-13(4) 






e-exp.[22] 


23.3 


2.12(12) 


2.42(30) 


-12(1) 






g-exp.[22j 


23.64(7) 


2.053(8) 


2.47(25) 


-25(18) 






<?-exp.[24] 


23.71 


2.060 










Approx. RG[25] 




1.938 


2.505 


-12.599 


10.902 




Approx. RG[26] 


20.72(1) 


2.063(5) 


2.47(5) 


-19(1) 






Approx. RG[27j 


28.9 


1.92 


2.17 








MC Ising sc[52] 


23.3(5) 


2.72(31) 










MC Ising sc[53] 


24.5(2) 


3.24(24) 











spin-s Ising models, we have taken advantage of the variable transformation eq. (l38p . to reduce 
the uncertainties of the estimates and the spread among the central values for different spins. We 
have always taken care that the uncertainties of our results allow for the errors of the critical 
temperatures listed in Table IVl and, whenever the variable transformation eq. (|38p is performed, 
also for the error of the exponent 9. 

In Figs. [6l...Il0]we have plotted vs the spin, our estimates of the quantities , , rt(j, rt2 
and for the Ising models of spin s = 1/2, ...3 on the sc and the bcc lattices. In these figures we 
have reported, in correspondence with the conventional value s = of the abscissa, also our results 
for the scalar model in the case of the sc lattice with quartic coupling g = 1.1 and, in the case of 
the bcc lattice, with coupling g = 1.85. The set of estimates shows good universality properties 
and moderate relative uncertainties which slowly grow with 2n. In the worst case, that of r^, the 
uncertainties are generally not larger than 15%. 

In Table HXl we have collected our final estimates of the ratios r^, rt,--- Tj*^, obtained either 
by simply choosing our result for the scalar field system on the bcc lattice, as in the case of the 
lowest-order ratios, or from a weighted average of the estimates on the sc and bcc lattices for the 
same system, as in the case of the largest-order ratios. Our values are compared with the estimates 
already obtained in the recent literature by various methods, including the analysis of significantly 
shorter HT expansions. 

V. CONCLUSIONS 

For a wide class of models in the 3D Ising universality class, we have described properties of the 
higher susceptibilities on the HT side of the critical point, which are relevant for the construction of 
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approximate representations of the critical ES. We have based on high-temperature and low- field 
bivariate expansions, that we have significantly extended or computed "ex novo". The models 
under scrutiny include the conventional Ising system with spin s = 1/2, the Ising model with spin 
s > 1/2 and the lattice scalar field, defined on the three-dimensional sc and bcc lattices. In this 
paper our HT data have been used to improve the accuracy and confirm the overall consistency 
of the current description of these models in critical conditions, by testing simple predictions of 
the scaling hypothesis as well as the validity of the universality property of the gap exponent and 
of appropriate combinations of critical amplitudes. Some of these tests are presently feasible only 
within a series approach. Our main result is a set of more accurate estimates of the first three 
already known parameters and a computation of two additional ratios, which enable us to 
formulate an update of the parametric form of the ES. 

At the order of expansion reached in our study, we still observe a small residual spread of the 
estimates of the gap exponent and of the ratios r^, around the predictions of asymptotic scaling 
and universality. This fact is readily explained by the obvious limitations of our numerical analysis: 
namely the still relatively moderate span of our expansions, in spite of their significant extension, 
the notoriously slower convergence of the expansions in the case of the sc lattice and the incomplete 
allowance of the non-analytic corrections to scaling by the current tools of series analysis. 
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FIG. 1: A plot vs the parameter p^, of the universal combinations of critical amplitudes 
C2^/C2" (upper curve), C^i?''~^/i?f (middle curve) and C^S^/(C^)^(lower curve) obtained from 
the truncated polynomial approximation of l{6) eq. (|24p . The computation is based on the coeffi- 
cients F2n-i with n = 1, ..7, estimated in this work. For convenience, the curves are normalized to 
their minimum values. 
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FIG. 2: The sequences of modified ratio approximants (MRAs) of the critical point for the scalar 
field model with g = 1.1 on the sc lattice, plotted vs l/r^"*"^. Here r is the order of the approximant 

and 6 is the exponent of the leading correction to scaling. We have normalized the MRAs to the 
estimated value of Kc- The MRAs are obtained from the HT expansions of X2{K) (stars), Xi{^) 
(circles), x&{K) (triangles), x?>{K) (rhombs), xio{K) (rotated squares), xi2{K) (squares), Xi4(i^) 
(double triangles), xi&{^) (rotated triangles), xis{K) (crossed circles), X2o{K) (crossed squares), 
X22{K) (crossed triangles). The symbols representing the MRAs are connected by straight lines as 
an aid to the eye. Small vertical segments on the third, fourth and fifth curve from above, indicate 
the order at which our extension of the x&{K), xsiK), and xio{K) series begins to contribute to 
the MRAs. The six lowest curves refer to higher susceptibilities for which no data exist in the 
literature. 
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FIG. 3: As a simple consequence of the scaling hypothesis, the exponent differences — '~i2n — 
72n-2 should not depend on n and equal 2A. Here they are obtained by forming second-order DAs 
of the ratios X2n{I^) / X2n-2{K) for n = 2,3,4, ...,11, obtained from the expansions of the hcc lattice 
Ising model with spin s = 1/2 (circles), s = 1 (triangles), s = 3/2 (rhombs), s = 2 (stars), s = 5/2 
(squares) and s = 3 (crossed circles). For each value of n, the various symbols have been slightly 
shifted apart only to avoid cluttering and keep the uncertainty of each estimate visible. The dashed 
horizontal line represents the estimated^ value 2A = 3.1276(8) of twice the gap exponent. The 
continuous horizontal lines indicate a deviation of 0.5% from the expected central value. 
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FIG. 4: Same as FigjSl In this case, the differences Dn = J2n — 72n-2 have been computed from 
the HT expansions of the higher susceptibihties for the lattice scalar field theory with values of 
the quartic coupling, g = 1.1 on the sc lattice (black squares) and g = 1.85 on the bcc lattice 
(black circles). These values of g are chosen to minimize the leading corrections to scaling. As 
in Fig 121 for each value of n, the symbols have been slightly shifted apart. The dashed horizontal 
line represents the estimated^ value 2A = 3.1276(8) of twice the gap exponent. The continuous 
horizontal lines indicate a deviation of 0.5% from the central value. 
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FIG. 5: The renormalized coupling constant ^4 for the Ising model with spin s on the bcc lattice 
(circles) and on the sc lattice (squares) vs the value s of the spin. The HT expansions of the 
Ising systems have been subjected to the variable transformation eq. (|38p . For comparison, we have 
also computed for the scalar model with the value g = 1.1 of the quartic coupling on the sc 
lattice (black square) and with g = 1.85 on the bcc lattice (black circle). The latter estimates are 
plotted with conventional abscissas near zero. In all cases, the symbols have been slightly shifted 
apart to avoid superpositions and to keep the uncertainty of each estimate visible. The dashed 
horizontal line represents the value g^ = 23.56(3) estimated in Ref.[r^. The continuous horizontal 
lines indicate a deviation of 0.2% from the central value. 
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FIG. 6: The ratio rg for the Ising model with spin s on the bcc lattice (circles), for the Ising model 
with spin s on the sc lattice (squares) vs the spin. The expansions for the Ising systems have been 
subjected to the variable transformation eq. (j38p . The symbols have been slightly shifted apart to 
avoid superpositions and to keep the uncertainties of each estimate visible. For comparison, we 
have also computed for the scalar model with g = 1.1 on the sc lattice (black square) and 
with g = 1.85 on the bcc lattice (black circle). The latter estimates are plotted with conventional 
abscissas near zero. 
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FIG. 7: Same as FigjGl but for r^. 
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FIG. 8: Same as FiglU but for r 
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FIG. 9: Same as FiglU but for r 
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FIG. 10: Same as FiglU but for r 
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